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0> ' Abstract 



We analyze the implications of the relation Trj^ = 0, which is customarily 

Q I assumed in practical lattice calculations. On the basis of the finite dimensional 

' representations of the Ginsparg- Wilson algebra, it is shown that this relation reflects 

• the species doubling in lattice theory; topological excitations associated with species 

CSJ , doublers, which have eigenvalue 2/a, contribute to Trj^ without any suppression. 

In this sense, the relation Tr'j^ = is valid only when we allow the presence of 

^ ' unphysical states in the Hilbert space; this statement is also valid in the Pauli-Villars 

CN . regularization. If one eliminates the contributions of the unpysical states, the trace 

I Tr'j^ is replaced by TrT^ = rr75(l — ^aD) which gives rise to the Pontryagin index, 

■ to be consistent with the continuum analysis. 

I ! 

as _ 

^ ■ 1 Introduction 



Recent developments in the treatment of fermions in lattice gauge theory led to a better 
understanding of chiral symmetry [1]- [9]. These developments are based on a lattice Dirac 
^ ■ operator D which satisfies the so-called Ginsparg- Wilson relation [1] 

K^: I5D + D-f5 = aD-i^D (1.1^ 

\ where 75 is a hermitian chiral Dirac matrix. An explicit example of the operator satisfying 
(1.1) and free of species doubling has been given by Neuberger[2]. The operator has also 
been discussed as a fixed point form of block transformations [3]. The relation (1.1) led 
to the interesting analyses of the notion of index in lattice gauge theory [4]- [9]. 
The index relation is often written as [4] [5] 

— Trj^-aD = — U- (1.2) 
where n± stand for the number of normalizable zero modes in 

DiPn = (1.3) 

with simultaneous eigenvalues j5<fn = ±V^n- 
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In the continuum path integral treatment of chiral anomaly, the relation 

Tr75 = n+ — n_ (1.4) 

in a suitably regularized sense (with help of a local version of the index theorem[10]) plays 
a fundamental role[ll]. On the other hand, it is expected that the relation 

rr75 = (1.5) 

holds on a finite lattice. As Chiu correctly pointed out [12], this relation (1.5) leads to an 
interesting constraint 

Tr75 = n+-n_ + A^+-A^_ = (1.6) 

where N± stand for the number of eigenstates D(pn — (2/a)</7„ with 75</?„ = ^fn, respec- 
tively. 

On the basis of the work [13], we here argue that Tr75 = implies the inevitable con- 
tribution from unphysical (would-be) species doublers in lattice theory or an unphysical 
bosonic spinor in Pauli-Villars regularization. In other words, Tr^^ — cannot hold in 
the physical Hilbcrt space consisting of physical states only, and the continuum limit of 
Tr75 = is not defined consistently, as is seen in (1.6). It is also shown that the failure 
of the decoupling of heavy fermions in the anomaly calculation is crucial to understand 
the consistency of the customary lattice calculation of anomaly where Tr^^ = is used. ( 
The continuum limit in this paper stands for the so-called "naive" continuum limit with 
a ^ 0, and the lattice size is gradually extended to infinity for any finite a in the process 
of taking the limit a 0.) 



2 Finite dimensional representations of the Ginsparg 
- Wilson algebra 

We start with the finite dimensional representations[13] of the basic algebraic relation 
(1.1). A construction of the operator D, which satisfies the Ginsparg- Wilson relation 
on a finite lattice, by using a corresponding operator D on an infinite lattice has been 
discussed in Ref.[14]. We first define an operator 

r5 = 75(l-^a^) (2.1) 
which is hermitian and satisfies the basic relation derived from (1.1) 

r575l^ + l^DT, = 0. (2.2) 
Namely, plays a role of 75 in continuum theory. This relation suggests that if 

75L>0„ = X„(f)n, (0n, 0n) = 1 (2.3) 

for the hermitian operator 75!^, then 

j5D{r,(Pn) = -A„(r5</.„). (2.4) 
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Namely, the eigenvalues A„ and — A„ are always paired if A„ 7^ and (Fs^^, Fs^^) 7^ 0. 
The inner product 0jj</)n = 0n) = Y.x '^'^'t>ni^)^n{x) is defined by summing over all the 
lattice points, which are not explicitly written in 
We evaluate the norm of Fs^^ 

ct ct 

(Fs^n, Fg^n) = (75 - -75!:)) (75 - -75i:))0„) 

= (1 - ^75(75/^ + i^75) + j(75/^)')0n) 

= j(75i^)')0n) 

= (1-|a„)(1 + |a„). (2.5) 



Namely 0„ is a "highest" state 



Fs^n = (75 - S75i^)</'n = (2.6) 



if (1 — f A„)(l + |An) = for the Euclidean 5'0(4)- invariant positive definite inner 
product {(i)n,<f>n)- Wc thus conclude that the states with A„ = ±^ are not paired by 
the operation F50„ and are the simultaneous eigenstates of 75, 750„ = ±0n respectively. 
These eigenvalues A„ are also the maximum or minimum of the possible eigenvalues of 
75D. This is based on the relation 

4>ilb(t>n = ^K<Pi<Pn = (2.7) 

for A„ 7^ 0, which is derived by sandwiching the relation(l.l) by 0]j75 and 0„. Namely, 

1^1 = \(pH<Pn\ < ||0n||||750n|| =1. 

For the vanishing eigenvalues 75-D0n = 0, one can show that '~^^D'~^<^(l)n = by using 
the relation (1.1). Thus one can choose the chiral eigenstates 75-D[(l ± 75)/2]0„ — 0. 
Namely 

lf,D(t)n = 0, 750„ = ±0„. (2.8) 
All the remaining normalizable states with < |A„| < 2/a, which appear pairwise 
with A„ = ±|A„| ( note that F5(F50„) = (1 — (aA„/2)^)0„ oc 0„ for |aA„/2| 7^ 1 ) , satisfy 
the relations 

0|^^50n = 0, 

/+ 1 ctA^ 

<t>Ll^<t>n = for Xm 7^ A„, A^A„ > 0. (2.9) 

These states 0„ cannot be the eigenstates of 75 as |aA„/2| < 1. 

The index (1.2) in the present context is most naturally defined as 

n 

An=0 0<|A„|<2/a X„=±2/a 

= E 475</'n = n+ - n_ (2.10) 

A„=0 



On the other hand, the index relation commonly written in the form [4] -[6] 

n 

= -fEAn = -A^H- + A^- (2.11) 

is saturated by the states iV±, where N± stand for the number of isolated (un-paired) 
states with A„ = ±2/a and 750„ = i^ru respectively. 

The relation Trj^ = 0, which is expected to be valid on a finite lattice, leads to ( by 
using (2.7)) 

n 

A„=0 A„^0 
A„=0 A„^0 ^ 

= ^ ^A„ = 0. (2.12) 

In the last line of this relation, all the states except for the states with A„ = ±2/a cancel 
pairwise for 7^ 0. We thus obtain a chirality sum rule n+ — n_ + iV+ — — [12] or, 

n+ + Ar+ = n_ + 7V_. (2.13) 

These relations show that the chirality asymmetry at vanishing eigenvalues is balanced 
by the chirality asymmetry at the largest eigenvalues with |A„| = 2/ a. 

Those properties we analyzed so far in this Section hold both for non-Abelian and 
Abelian gauge theories. We did not specify precise boundary conditions, since our anal- 
ysis is valid once non-trivial zero modes appear for a given boundary condition. For an 
Abelian theory, one needs to introduce the gauge field configuration with suitable bound- 
ary conditions, which carries a non-vanishing magnetic flux , to generate a non-trivial 
index — n_ [14]. Our analysis of the index in this Appendix is formal, since it is 
well known that the Ginsparg- Wilson relation (1.1) by itself does not uniquely specify 
the index or the coefficient of chiral anomaly for a given gauge field configuration [15] .A 
suitable choice of parameters in the overlap operator [2], for example, uniquely specifies 
the index. 

To summarize the analyses of the present Section, all the normalizable eigenstates 4>n 
of the hermitian 75 D on a finite lattice are categorized into the following 3 classes: 
(i) n± states. 



75L>(/)„ = 0, 75(/)„ = ±(/)n, (2.14) 

(n) N± states(r50„ = 0), 

2 

^^D(j)n = ±-0„, 750„ = ±0„, respectively, (2-15) 
a 
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(iii) Remaining states with < |A„| < 2/a, 

75^0n = K4>n, 75^(r50n) = -A„(r50„), (2.16) 

and the sum rule n+ + iV+ = n_ + holds. 

All the n± and N± states are the eigenstates of D, D(j)n = and -D0„ = (2/a)0„, 
respectively. If one denotes the number of states in (iii) by 2No, the total number of 
states (dimension of the representation) N is given hy N = 2(n+ + + Nq), which is 
expected to be a constant independent of background gauge field configurations. 



3 The nature of the N± states with T^cpn = 

In the previous section we have seen that the consistency of the relation Tr'-fr, = requires 
the presence of the N± states for an operator 75D satisfying (1.1) on a finite lattice. 
We want to show that the N± states are topological excitations associated with species 
doublers. For this purpose, we start with the conventional Wilson operator Dw 

Dw{n, m) = iYCf^{n, m) + B{n, m) - ^moSn,m, 
C^{n,m) = ^[Sm+ij.,nUij,{m) - S^^n+t,Ul{n)], 

B(n, m) = — Y^[2Sn,m - 5m+tji,nU ,j.{m) - 5m,n+,xUl{n)], 
za ^ 

U^{m) = exp[iagA^{m)], (3.1) 

where we added a constant mass term to Dw. Our matrix convention is that 7'* are 
anti-hermitian, (7'')^ = —J^, and thus = j^CulnjUi) is hermitian 

(3.2) 

In the case of a very explicit construction of species doublers is known. For a square 
lattice one can explicitly show that the simplest lattice fermion action 

s^i)ig:i) (3.3) 

is invariant under the transformation [16] 

V'' = TV', i)' = ^r-^ (3.4) 

where T stands for any one of the following 16 operators 

1, T{r2, T,Ts, T,T^, TsTa, T^T^, T^T^, T^T^T^T^, (3.5) 

and 

Ti, Ta, Tg, T4, TiTaTa, 72^3^4, T^T^^u T^T^T^. (3.6) 
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The operators are defined by 

7; = 7^75exp(i7rxVa) (3.7) 

and satisfy the relation 

T^T, + T,T^^2S^, (3.8) 

with = = T^^ for anti-hermitian 7^. We denote the 16 operators by 7^, n = 0~15, 
in the following with % — 1. By recalling that the operator 7)t adds the momentum tt/o 
to the fermion momentum /c^, we cover the entire Brillouin zone 

_ ^ < A;^ < ^ (3.9) 

by the operation (3.4) starting with the free fermion defined in 

-^<.,<|^. (3.10) 

The operators in (3.5) commute with 75, whereas those in (3.6) anti-commute with 75 and 
thus change the sign of chiral charge, reproducing the 15 species doublers with correct 
chiral charge assignment; Z]^lo(~l)"'75 — 0- 

In a smooth continuum limit, the operaor produces p for each species doubler with 
alternating chiral charge. The relation Tr'-f^ = for the operator ^ is consistent for 
any background gauge field because of the presence of these species doublers, which are 
degenerate with the physical species in the present case. 

The consistency of Tr'-f^ = is generally analyzed by means of topological properties 
and thus it is best described in the nearly continuum limit. To be more precise, one may 
define the near continuum configurations by the momentum carried by the fermion 

- —e<k.< —e (3.11) 
2a ^ ^ - 2a ^ ' 

for sufficiently small a and e combined with the operation 7^ in (3.5) and (3.6). To identify 
each species doubler clearly in the near continuum configurations, we also keep r/a and 
vriQja finite for a small [16], and the gauge fields are assumed to be sufficiently smooth. 
For these configurations, we can approximate the operator by 

= i ^ + M„ + O(e^) + O(a^A^) (3.12) 

for each species doubler, where the mass parameters M„ stand for Mq — and one of 

2r mo , . mo , . 

, (4,-1); , (6,1) 

a a a a 

(4.-1); (1.1) (3.13) 

a a a a 

for n = 1 ~ 15[2]. Here we denoted ( multiplicity, chiral charge ) in the bracket for species 
doublers. In (3.12) we used the relation valid for the configurations (3.11), for example, 

Dw{k) = 2^7 - + cos ak I,) 

= Yk,{l + 0{e')) + -0{e')-'^ (3.14) 

a a 
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in the momentum representation with vanishing gauge field. 

In these near continuum configurations, the topological properties are specified by the 
operator ]/) in D^/. The overlap operator D introduced by Neuberger[2] , which satisfies 
the relation (1.1), has an explicit expression 

aD^l- = 1 + Dw I ] (3.15) 

Vi^2 ^D^Dw 

where Dw — is the Wilson operator. For the near continuum configurations spec- 

ified above in (3.11), one can approximate 

D = f;(i/a)[l + (i^ + M„)^=L==]|n)(n|, 

n=0 ^P' + M^ 

= ^(-l)"75(l/a)[l + (^^ + M0 , / ,J \^){M, 
n=o y/W+K 

75 = E(-l)"75|n)(n|. (3.16) 

n=0 

Here we explicitly write the projection \n){n\ for each species doubler. The operators in 
(3.16) preserve the Ginsparg- Wilson relation (1.1). 
The above expression of D in (3.16) shows that 

D(Pi = 0, 

Dcj>i = (3.17) 
a 

for the physical species and the unphysical species doublers, respectively, if one uses the 
zero-modes in 

P(f)l = Xn(f>n- (3.18) 

Note that Mq < and the rest of M„ > in (3.13) and (3.16) [2]. We also note that 
(pi can be a simultaneous eigenstate of 75 only for p(t)i = 0. Namely, the N± states with 
eigenvalue 2/a of D in fact correspond to topological excitations associated with species 
doublers. This means that the multiplicities of these N± are quite high due to the 15 
species doublers, although they satisfy the sum rule + = n_ + N_: This sum rule 
itself is a direct consequence of (3.17), (3.18) and Z^„(— 1)^75 = if one recalls ^1750; = 
for \i ^ 0. The relation (3.17) with zero-modes in (3.18) shows that the index relation 
(2.10) stands for a lattice version of the Atiyah-Singer index theorem. 

We here briefiy sketch an evaluation of the anomaly, since it justifies our analysis based 
on the effective expressions in (3.16). For an operator 0{x,y) defined on the lattice, one 
may define 

Omn = J2<l>*mi^)0{x,y)My), (3-19) 
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and the trace 



TrO = Y^Onn 

n 

= EEC(^)O(x,y)0„(y) 

n x,y 

= E(EC(^)O(a:,?/)0„(z/)). (3.20) 



X n,y 



The local version of the trace (or anomaly) is then defined by trO{x,x) = 
J2n,y (l^n{^)0{x,y)(l)n{y)- For the operator of our interest, we have 

trj,{l - ^D){x) = E tr ^^e-^'^T^-'j,^^=V^ (3.21) 
2 2 7-^ (27r)4 \Jd\^Dw 

where we used the plane wave basis defined in (3.10) combined with the operation 7^. 
We also used the chiral charge assignment (3.16) and a short hand notation Oe^^^ — 

We first take the a — >■ limit of (3.21) with all M„,n = ~ 15, kept fixed and then 
take the limit |M„| — > oo later. For fixed M„ ( to be precise, for fixed mo/a and r/a ), 
one can confirm that the above integral (3.21) for the domain < |A;^| < ^ vanishes 
( at least ) hnearly in a for a —> 0, if one takes into account the trace with 75. See also 
Refs.[7][9]. In the remaining integral 

- ^ U-mr /: ^,e-^'^j,T-^D^^^%e^'^ (3.22) 
2 „=o -^-^a^ (27r)^ ^JDlDw 

one may take the limit a — > ( and — > 00 ) with letting e arbitrarily small. By taking 
(3.12) into account, one thus recovers the local anomaly [8] 

2„^; ^ y-»(27r)^ ^ 

1 /-oo 1 15 -1 

2 7-^ (27r)4 ''^^^2/^2 + 1 ^jp^/Ml-rl 

= ^ire'^'^^^i^M.i^a/? (3.23) 

for \Mn\ 00. The domain in (3.11) with arbitrarily small but finite e thus correctly 
describes the topological aspects of the continuum limit in the present prescription. 

Here we went through some details of the anomaly calculation to show that the in- 
terpretation of the N± states in (2.15) as topological excitations associated with species 
doublers, as is shown in (3.17), is also consistent with the local anomaly calculation. As 
for a general analysis of chiral anomaly in the overlap operator, see Ref.[17]. 
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3.1 General lattice Dirac operator and Tr^y^ = 



We expect that our analysis of Tr^^ = 0, namely its consistency is ensured only by the 
presence of the would-be species doublers in the Hilbert space, works for a general lattice 
Dirac operator, since any lattice operator contains as an essential part. For the smooth 
near continuum configurations, the lowest dimensional operator is expected to specify 
the topological properties. 

We also note that the Pauli-Villars regularization in continuum theory can be analyzed 
in a similar way. The Pauli-Villars regulator is defined in the path integral by introducing 
a bosonic spinor into the action 

S = J d'^x[i/j{i p - m)tl) + M)0]. (3.24) 

The Jacobian for the global chiral transformation then gives rise to the graded trace[ll] 

Tr75 = Tr^75 - Tr^j^ = 0. (3.25) 

The relation Tr^r^ = is thus consistent with any topologically non-trivial background 
gauge field because of the presence of the unphysical regulator 4>. This 4> is analogous to 
the species doublers in lattice regularization. 



4 Implications of the present analysis 

We have shown that the consistency of Tr^r, = for topologically non-trivial background 
gauge fields requires the presence of some unphysical states in the Hilbert space. Coming 
back to the original lattice theory defined by 

S = jj^Dijj (4.1) 

with D satisfying the relation (1.1), one obtains twice of (2.10) as a Jacobian factor for 
the global chiral transformation [5] 6t/j — ^675(1 — |-D)V' and St/j — ■0ie(l — |-D)75, which 
leaves the action (4.1) invariant. One can rewrite (2.10) as 

Tr75(l - |d) = Tr75(l - |d) = fn^ = n+ - n. (4.2) 

where the modified trace Tr is defined by truncating the unphysical N± states with 
A„ = ±2/a . Without the N± states, Tvy^^D — since the eigenvalues A„ of ^^D with 

An 7^ 0, ±2/a appear always pairwise at ±|A„|. 

If one takes a smooth continuum limit of Tr'-f^ = — n_ in (4.2) , one recovers the 
result of the continuum path integral (1.4). If one considers that Tr75 is too abstract, 
one may define it more concretely by 

T.7.(l-|i^)/(^) = T.7a(l-|i.)/(^) 

= Tr75/(^^) = n+-n_ (4.3) 
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for any f{x) which rapidly goes to for x — > oo with /(O) = 1. This relation suggests 
that we can extract the local index [10] ( or anomaly) by 

tn.i^-^D)fi^^^)i^) (4.4) 

which is shown to be independent of the choice of f{x) in the limit a and leads to 
(3.23) by using only the general properties of D [8]. We thus naturally recover the result 
of the continuum path integral[ll]. As for an interesting algebraic analysis of the local 
anomaly tr'y^{l — f-D), see [18]. 

As for a more practical implication of Tvf^ = in lattice theory, one may say that 
any result which depends explicitly on the presence of the states N± is unphysical. It is 
thus necessary to define the scalar density ( or mass term ) and pseudo-scalar density in 
the theory (4.1) by [19] [6] 

S{x) = IpLlpR + ijRlpL = IPi^ - ^^)^' 

P{x) = ^j^L^R-^i^RJjJL^^l5(l-^D)^|J. (4.5) 
Here we defined two independent projection operators 

P± = ^(1±75) 

P± = ^(1±75) (4.6) 

with 75 = 75(1 — aD) which satisfies 7I = 1 [6]. The left- and right- components are then 
defined by 

^L,R^'^P±, i^R,L^P±i^ (4.7) 

which is based on the decomposition 

D = P+DP_ + P_DP+. (4.8) 

The physical operators S{x) and P{x) in (4.5) do not contain the contribution from the 
unphysical states N± in (2.15). In the spirit of this construction, the definition of the 
index by (4.3) which is independent of unphysical states N± is natural. In particular, all 
the unphysical species doublers ( not only the topological ones at 2/a ) decouple from the 
anomaly defined by (4.4) in the limit a — > with fixed M. 

The customary calculation of the index ( and also anomaly ) by using the relation [4]- 
[7] [9] 

rr75(l - ^D) = Tr{-^-f5D) ^ n+ - (4.9) 

by itself is of course consistent, since one simply includes the unphysical states N± in 
evaluating Tr^^ — 0, and consequently one obtains the index T'r(— |75L>) from the un- 
physical states N± only. We after all know that the left-hand side of (4.9) is independent 
of A^±. 
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Rather, the major message of our analysis is that the continuum hmit of Tr'~)^ = 
in (1.6) ( unhke the relation Tr75 = n+ — n_ ) cannot be defined in a consistent way 
when the (would-be) species doublers disappear from the Hilbert space. It is clear from 
the expression of Trj^ = in (1.6) that the a —>■ limit of Trjr, = is not defined 
consistently. One may then ask how the calculation of local anomaly on the basis of (4.9) 
could be consistent in the limit a — > if Tr75 = is inconsistent? A key to resolve this 
apparent paradox is the failure of the decoupling of heavy fermions in the evaluation of 
anomaly. The massive unphysical species doublers do not decouple from the anomaly , 
as is seen in (3.23), for example. If one insists on Tr'f^ = in the continuum limit, one 
is also insisting on the failure of the decoupling of these infinitely massive particles from 
Tr75 = 0. The contributions of these heavy fermions to the anomaly and to Tr'y^ = 
precisely cancel, just as in the case of the evaluation of global index in (4.9). Namely, 
the local anomaly itself is independent of these massive species doublers in the continuum 
limit, as is clear in (4.4). It is an advantage of the finite lattice formulation that we can 
now clearly illustrate this subtle cancellation of the contributions of those ultra-heavy 
A^-t-states to Tr'-f^ = and anomaly on the basis of (1.6). 

The unphysical nature of the A^±-states has been recently clarified from a different 
view point in Re. [20] 
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